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)Q Abstract 

o ■ 

""^5 ' We studied extended X-ray emission from the Carina Nebula taken with the Suzaku CCD camera XIS 

on 2005 Aug. 29. The X-ray morphology, plasma temperature and absorption to the plasma are consistent 
with the earlier Einstein results. The Suzaku spectra newly revealed emission lines from various species 
including oxygen, but not from nitrogen. This result restricts the N/O ratio significantly low, compared 
with evolved massive stellar winds, suggesting that the diffuse emission is originated in an old supernova 
remnant or a super shell produced by multiple supernova remnants. The X-ray spectra from the north and 
south of 77 Car showed distinct differences between 0.3—2 kcV. The south spectrum shows strong L-shell 
, ' lines of iron ions and K-shcll lines of silicon ions, while the north spectrum shows them weak in intensity. 

, This means that silicon and iron abundances are a factor of 2—4 higher in the south region than in the 

' north region. The abundance variation may be produced by an SNR ejecta, or relate to the dust formation 

around the star forming core. 
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1. Introduction 

Not a few giant Hll regions emit extended soft X-rays 
with fcT~0.1-0.8 keV, log Lx ^33-35 ergs s _1 , and extent 
~1-10 3 pc (RCW 38: Wolk et al. 2002, M17: Townsley 
et al. 2003, NGC 6334: Ezoe et al. 2006, Carina Nebula: 
Seward et al. 1979; Seward & Chlebowski 1982; Evans 
et al. 2003, extragalactic Hll regions: Kuntz et al. 2003, 
Strickland et al. 2004). This high surface brightness can- 
not be accounted for by an extrapolation of the discrete 
source luminosity function, and thus is thought to arise 
from diffuse plasma. The required plasma temperature 
and thermal energy can be produced by collisions or ter- 
mination of fast winds from main-sequence or embed- 
ded young O stars (Seward & Chlebowski 1982; Townsley 



et al. 2003; Ezoe et al. 2006), but the extended emission 
is often observed from regions apart from massive stellar 
clusters. Another origin, such as an unrecognized super- 
nova remnant (SNR), cannot be ruled out. 

In principle, the origin of the diffuse emission can be 
determined by measuring its composition. For example, 
the plasma should be overabundant in nitrogen and neon 
if it originates from winds from a nitrogen-rich Wolf-Rayet 
star (WN), while it would be overabundant in oxygen if 
it arises from a Type II SNR. The temperature of the 
plasma, typically a few million degrees, make soft X-ray 
band studies highly desirable, because of the presence in 
this band of strong lines from these elements, plus carbon, 
silicon and iron. 

The Carina Nebula, which contains several evolved and 
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main-sequence massive stars such as rj Car, WR 25 and 
massive stellar clusters such as Trumpler 14 (Tr 14), emits 
soft diffuse X-rays 10-100 times stronger than any other 
Galactic giant Hll region (Lx ~10 35 ergs s _1 , Seward 
et al. 1979). The high surface brightness made possi- 
ble the discovery of the diffuse emission by the Einstein 
Observatory in the late 1970's. The Einstein observations 
revealed that the diffuse emission tends to be associated 
with optically bright regions containing massive stars, and 
is strongest in regions bordering the western part of the 
V-shaped dust lane (Seward & Chlcbowski 1982). Recent 
Chandra observations provided a point source free mea- 
surement of the diffuse flux (Evans et al. 2003), and sug- 
gested presence of a north-south Fe and Ne abundance 
gradient (Townsley 2006). 

No detailed abundance measurement of the diffuse emis- 
sion from the Carina Nebula has been made. The X-ray 
CCD cameras (XISs: X-ray Imaging Spectrometer) on- 
board the Suzaku observatory have the best spectral reso- 
lution for extended soft X-ray emission. It is thus suitable 
for measuring the abundance of key elements in the diffuse 
plasma. To compensate for the limited angular resolution 
of Suzaku, we have also utilized the archival XMM-Newton 
(Jansen et al. 2001) and Chandra (Weisskopf et al. 2002) 
data from the Carina Nebula to investigate the spatial 
distribution of the diffuse emission and the contamination 
from X-ray point sources. 

2. Observations &; Data Reduction 

Suzaku (a.k.a. Astro-E2, Mitsuda et al. 2006) is 
the 5th Japanese X-ray observatory, developed by a 
Japanese-US collaboration, and launched on 2005 July 
10. The observatory is equipped with five thin-foil X-Ray 
Telescopes (XRT) with half power diameter (HPD) of ~2' 
(Serlemitsos et al. 2006). At their focal planes are four 
CCD focal plane detectors (XIS, Koyama et al. 2006b) 
and an X-ray calorimeter (XRS, Kelly et al. 2006), which 
unfortunately malfunctioned a month after the launch. 
The observatory also has a hard X-ray detector (HXD, 
Takahashi et al. 2006; Kokubun et al. 2006). The XIS 
comprises four CCD cameras, XIS0 through XIS3, three 
of which (XIS 0, 2 and 3) use front-illuminated (FI) CCD 
chips while XIS1 utilizes the back- illuminated (BI) tech- 
nology. Compared with CCD detectors on earlier X-ray 
observatories, both the BI and FI chips have better spec- 
tral resolution at energies below ~1 keV: their FWHM 
energy resolution at 0.5 keV is ~40 eV (FI) and ^50 eV 
(BI), with negligible low-energy tails. 

Suzaku observed the Carina Nebula twice, on 2005 Aug. 
29 and 2006 Feb. 3, with similar pointings, putting -q Car 
at the center of the XIS FOV (Figure 1). In both observa- 
tions, the XISs were operated with the Normal mode. The 
first observation was 16 days after the XIS first light (2005 
Aug 13) when contamination on the optical blocking filter 
was almost negligible, while the second observation, 174 
days after the first light, suffered significant degradation 
of the soft response by progressive contamination on the 
XIS of a factor of ~2 between 0.4-0.5 keV on axis. We 



thus did not analyze the second observation in this paper. 
We also did not use the HXD data, whose X-ray signals 
are mostly from rj Car. 

We used the HEAsoft 1 analysis package ver. 6.0.6 for 
data analysis. For the version 0.7 2 unscreened data, we 
screened out hot and flickering pixels with cleansis im- 
plemented in the xselect and removed exposure frames 
with telemetry saturation. We excluded the data taken 
during the South Atlantic Anomaly (SAA) passage and 
within 256 sec after the passage. We selected data taken 
while the geomagnetic cut-off rigidity (COR) was larger 
than 4 GV, elevation from above the night-time Earth 
horizon (ELV) was larger than 5°, and elevation above 
the sun-lit Earth rim (DaY Earth ELeVation: DYE_ELV) 
was larger than 10° for the FI data and 20° for the BI. The 
tight DYE_ELV filtering criterion for the BI was intended 
to minimize contamination by neutral nitrogen and oxy- 
gen lines from the sunlit terrestrial atmosphere. There 
were no remarkable particle background flares that de- 
graded the data quality. The net exposure of the each 
FI chip and the BI were ~60 ksec and 49.0 ksec, respec- 
tively. In spectral fits, we used the XIS response matri- 
ces released on 2006-02-13 and generated auxiliary files 
with xissimarf gen-2006-05-28, which considers addi- 
tional absorption due to the XIS contamination. 

3. Spatial Distribution of the X-ray Sources 

To understand the spatial distribution of the diffuse X- 
ray emission and X-ray point sources in the field, we sur- 
veyed archival XMM-Newton data-sets aimed at r\ Car. 
There are in total 11 observations performed between 2000 
and 2003 3 . Figure 1 shows a true color image accumu- 
lated from all the MOS data taken with the prime full 
window mode (red: 0.2—0.6 keV, green: 0.6—1.2 keV, 
blue: 1.2—10 keV). The image depicts several bright X- 
ray point sources: n Car (LBV star), WR25 (Wolf-Rayet 
star), HD 93250 and HD 93205 (03 stars), and Tr 14 
(massive stellar cluster), as well as multiple weak point 
sources especially concentrated between r\ Car and WR 25. 
The image also clearly shows two kinds of apparently ex- 
tended emission, seen as red and green. The "red" emis- 
sion — relatively strong below 0.6 keV — is located be- 
tween Tr 14, WR 25 and 77 Car. As is suggested by Seward 
& Chlebowski (1982), the emission is stronger at the west 
side with a sharp cut-off at the boundary of the V-shaped 
optical dark lane. The "green" emission — relatively 

1 http:/ /heasarc. gsfc.nasa.gov/docs/softwarc/lhcasoft/ 

2 The version 0.x process uses the earliest detector calibration and 
Suzaku software still under development, and was released only 
to the Suzaku team members (Science Working Group: SWG) 
to verify the data quality. The process does not perform some 
detailed calibration of the data, such as correction of the satellite 
wobbling and event time tagging optimized for the XIS opera- 
tional mode. The latest processing as of 2006 July is version 0.7, 
which achieved an absolute energy scale within ±5 eV below 1 
keV and ±0.2% at the iron Kq energy for the XIS data. 

3 Sequence ID: 112580601, 112580701, 145740101, 145740201, 
145740301, 145740401, 145740501, 160160101, 160160901, 
145780101, 160560101, 160560201, 160560301. 
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Fig. 1. Top left: XMM-Newton EPIC MOS image color-coded to represent the soft band (0.2-0.6 keV) to red, medium band 
(0.6-1.2 keV) to green and hard band (1.2-10 keV) to blue. The bar-dot lines show the XIS FOV in the 2005 Aug. 29 observation. 
Top right: Suzaku BI image in the 2005 Aug. 29 observation. The solid white lines show event extraction regions, and the solid 
green lines show the Chandra FOV in the 1999 Sep. 6 observation. Bottom left: Integrated intensity map of the 12 CO emission 
(sky blue contour) on the optical image taken from the Digitized Sky Survey (gray scale image), which is a part of Figure 1 in 
Yonekura et al. (2005). The red spots show the bright X-ray sources introduced in the top left image. The yellow circle shows the 
XMM-Newton EPIC MOS FOV. Bottom right: Chandra image in the 1999 Sep. 6 observation. All the X-ray images are drawn with 
logarithmic scale and do not correct the vignetting effect. 
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Fig. 2. BI (top) and FI (bottom) spectra of the north (left) and south (right) regions. The black shows the source spectra and 
the grey shows the detector internal background, obtained from night Earth observations. The above labels demonstrate energies of 
emission lines detected or concerned with this result. 



strong above 0.6 keV — is located to the south of r\ Car, 
extending cast-west, where another optical dark lane and 
CO gas condensation are located (Yonckura et al. 2005, 
also see Figure 1). 

The Suzaku XIS FOV (Figure 1) covers most of the 
X-ray point and extended sources detected in the XMM- 
Newton image. Among the point sources, the XIS im- 
age resolves four bright X-ray sources, r/ Car, WR 25, 
HD 93250 and HD 93205. To minimize contamination, 
we excluded regions around these sources, within 4' from 
r\ Car and 2' from the other three sources. We then defined 
the "north" and "south" regions as in Figure 1 , to extract 
spectra of the "red" and "green" emission, respectively. 
We excluded an area between the north and south regions 
containing multiple weak point sources. Bright spots on 
the BI image at the upper and lower right corners are 
55 Fe calibration sources. We included these regions in the 
BI data to increase the photon statistics for the diffuse 
emission, since the calibration sources do not significantly 
contaminate the spectrum outside 5.9 and 6.4 keV, while 
we removed them from the FI data so we could investigate 



the Fe K line profile around 6—7 keV. 

4. Characteristics of the North and South 
Spectra 

We extracted all the XISO-3 spectra from both the north 
and south regions and combined the spectra of the FI chips 
(XIS0, XIS2 and XIS3), whose spectral responses are al- 
most identical. (Hereafter, we call the combined spectrum 
as "FI spectrum" .) We overlaid these FI and BI spectra of 
the north and south regions on spectra accumulated from 
observations of the night side of the Earth (Figure 2). 
Events detected during the night Earth observations orig- 
inate from particle background, X-ray fluorescence of Ni, 
Au and Si inside the camera body, 55 Fe calibration source 
and detector electronic noise. Accurate reproduction of a 
background spectrum requires proper accounting for the 
range of COR values encountered during the observation. 
For each of the two regions, we accumulated in detector co- 
ordinates for each field night Earth spectra corresponding 
to various COR values and constructed a spectrum from a 
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weighted average of these. The night Earth spectra match 
with the observed spectra at >10 keV for the FI with less 
than ^7.5% discrepancy, and >8 keV for the BI with less 
than ~2.5% discrepancy. A small discrepancy by ~10% 
is seen in lines from the 55 Fe calibration source, which de- 
cayed from the Carina Nebula observation on 2006 Aug 29 
to the span between 2005 September - 2006 May when the 
night Earth spectra were accumulated. This means that 
the internal background during this observation is well re- 
produced with the night Earth spectra, and we therefore 
use the night Earth data as a background in subsequent 
analysis. Both the north and south spectra show sub- 
stantial excess X-ray emission above the background level 
between 0.2—10 keV, including prominent emission lines. 

The north and south spectra show similar spectral fea- 
tures above 2 keV and below 0.3 keV. Above 2 keV, both 
spectra show emission lines from S XV, S XVI, Ar XVII and 
FeXXV ions, emitted from hot plasma with kT >2 keV, as 
well as a marginal detection of fluorescence from cold iron 
at 6.4 keV. The south region has ~60— 70% of the number 
of events from the north region, which just corresponds to 
the ratio of the effective area, and does not strongly de- 
pend on the photon energy. This means that the flux level 
and continuum slope are also similar between the north 
and south spectra above 2 keV. This similarity is seen for 
emission below 0.3 keV as well, though it is only measur- 
able with the BI. The north and south spectra show no 
remarkable variation below 0.3 keV and above 2 keV. 

In contrast, the spectra between 0.3—2 keV are con- 
spicuously different. Figure 3 shows an overlay of the 
BI spectra between 0.3-2 keV, with the north spectrum 
normalized by 64% to compensate for the effective area 
difference at 2 keV. A strong difference is seen between 
0.7 keV and 1.2 keV, which apparently is the source of 
the two colors of diffuse emission in Figure 1. The band 
in which the difference is found is dominated by emis- 
sion lines from the iron L-shell complex. Additionally, the 
south spectrum shows a stronger SiXIII line and possibly 
a stronger CVI line at ^0.39 keV, while the north spec- 
trum shows a clear excess at and below the O VII line. The 
south spectrum shows no clear evidence of lines from Ne IX 
or Ne X ions, but these could be hidden by the Fc L-shcll 
complex. Except for the energy bands with these emission 
lines, the spectral shapes look similar. This suggests that 
the differences represent an elemental abundance varia- 
tion, and not a temperature difference. We should note 
that though a small marginal peak at 0.5 keV looks like 
an emission line of hydrogen-like nitrogen, it is more likely 
produced by a combination of interstellar absorption and 
the instrumental oxygen edge at 0.54 keV. 

These results suggest that the spectra from both re- 
gions arc roughly divided into three components: emis- 
sion above 2 keV (hard component, hereafter HC) and 
below 0.3 keV (soft component, SC), which do not differ 
significantly between the north and south regions; and 
emission between 0.3 and 2 keV (medium component, 
MC) with remarkable variation in some emission lines. 
We see strong variation in emission lines from the Fe L- 
shell complex in the MC, but not even slightly in emission 
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Fig. 3. North (black) and south (red) spectra between 0.3—2 
keV overlaid. The north spectrum is normalized by 64%, to 
adjust the effective area of the south region at 2 keV. The 
above labels show detected emission lines, among which blue 
elements exhibit distinct differences between the north and 
south spectra. 

lines from FeXXV ions in the HC. This supports the idea 
that each component has a different origin. The most 
plausible origin of the SC is foreground emission; for the 
HC it is a combination of point sources, Cosmic X-ray 
Background (CXB) radiation, and Galactic Ridge X-ray 
Emission (GRXE); and for the MC, truly diffuse emission 
from the Carina Nebula. In the next section, we estimate 
these contributions in detail. 

5. Contamination from X-ray Sources 

Both spectra include emission from numerous point 
sources that are not resolved with Suzaku (see the XMM- 
Newton image in Figure 1), CXB, GRXE and Local Hot 
Bubble (LHB) emission. In this section, we estimate the 
contribution of these contaminants in the Suzaku spec- 
tra. We plot their estimated contribution in Figure 4, 
using point or extended source responses generated with 
xissimarf gen-2006-05-28. 

5.1. Point Sources 

In a Chandra observation in 1999, Evans ct al. (2003) 
detected 25 and 12 X-ray point sources from the north 
and south regions, respectively. They estimated the X-ray 
flux of each detected source from its count rate, assuming 
a typical stellar spectrum and extinction to the Carina 
Nebula (log T = 6.65 K, N R = 3xl0 21 cm" 2 ). However, 
several point sources in the field are found to have signif- 
icantly higher plasma temperature than the typical stel- 
lar coronal source (Albacete Colombo et al. 2003; Evans 
et al. 2004). This means that representing the contri- 
bution of unresolved point sources in the Suzaku spectra 
using a typical stellar spectrum scaled to the total point 
source flux calculated from tables in Evans et al. (2003) 
would be incorrect. We therefore analyzed the archival 
Chandra observation of the Carina Nebula. 

The Chandra observation we analyzed was performed 
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Fig. 4. Background subtracted BI spectra (left: north, right: south) with possible spectra of contamination sources. The point 
source spectra are scaled to the area coverage of the Suzaku source regions (17% to the north and 59% to the south). 



with the ACIS-I detector on 1999 Sep. 6 for 21.5 ksec (se- 
quence ID: 50, 1249), and has about twice the exposure 
of the data set used in Evans et al. (2003). In this ob- 
servation, the ACIS-I covered most of the XIS FOV (see 
Figure 1). We used the CIAO ver. 3.3 for the analysis. 
Using wavdetect on the 0.5—8 keV image and setting the 
detection threshold to a significance 10~ 6 , we detected 
49 and 25 sources in the north and south regions respec- 
tively. We extracted X-ray events of these point sources 
and background from a nearby source free region. Figure 5 
depicts the composite spectra of the point sources in the 
north and south regions. We generated response matri- 
ces and auxiliary files with mkrmf and mkwarf . The north 
spectrum shows weak neon and magnesium features. In 
contrast, the south spectrum appears featureless, but this 
could be due to low spectral resolution of this data taken 
in an early observing phase without CTI correction. Each 
spectrum was fitted using an absorbed 2-tcmpcrature (2T) 
model (WABS, Morrison & McCammon 1983; the APEC 4 
code, Table 1). A model of the south spectrum is repro- 
duced at above 90% confidence, while a model of the north 
spectrum is reproduced at slightly less than 90% confi- 
dence, possibly because of marginal line features around 
1 keV. Both best-fit models have high kT ^4 keV and 
consistent Nn values with the interstellar absorption to 
the Carina Nebula. The extremely low abundances are 
artifacts of an incompletely calibrated response for the 
early Chandra data and/or combining spectra from mul- 
tiple sources. 

The ACIS-I FOV does not cover -15% of the north re- 
gion and —37% of the south region. We therefore scaled 
the normalization of the spectra shown in Figure 4 to 
account for the partial coverage, though in the XMM- 
Newton image only a few point sources are detected in the 
uncovered portion of the north missing region and none 
in the south. The estimate still has an uncertainty in that 
variability of point sources must be taken into account, 
though we do not expect sources to vary in concert. The 

4 http:/ /cxc. harvard.edu/atomdb/sources_apec. html 



point source contribution to the Suzaku to the north and 
south regions inferred from our analysis of the Chandra 
data would be <l/3 of the total flux. Townsley (2006) 
suggests that the number of faint X-ray sources around 
the Tr 14 cluster could exceed one thousand; it is there- 
fore possible that the contribution from faint point sources 
is underestimated. 

5.2. CXB, GRXE, LHB and 77 Car 

For a CXB spectrum we use model Idl in Tabic 2 of 
Miyaji et al. (1998), which suggests a CXB flux fluc- 
tuation between regions of <30%. We assumed at 
1.3xl0 22 cm~ 2 , the Galactic Hi column density to the 
Carina Nebula (Dickey & Lockman 1990). The CXB flux 
contribution to the north and south spectra is <10%. 

To estimate the GRXE flux around the Carina Nebula 
(I, b) ~(287.°60, -0.°63) we referred to a 3-20 keV 
Galactic plane map taken with Rossi X-ray Timing 
Explorer (RXTE) (see the top panel of Figure 7 in 
Revnivtsev et al. 2006). In this map, the Carina Nebula 
is heavily contaminated by emission from ?y Car because of 
the limited RXTE spatial resolution. We therefore mea- 
sured X-ray flux at slightly higher (I —290°) and lower 
(I ~285°) Galactic longitudes along the Galactic plane, 
where bright X-ray point sources are apparently absent, 
and interpolated those values to infer a GRXE flux at 
the Carina Nebula of ~ 1.4xl0 _u ergs cm" 2 s _1 deg -2 
(3—20 keV). We assumed the GRXE spectral shape mea- 
sured by Ebisawa et al. (2005) at (I, b) ~(28.°5, 0.°0) and 
reduced the normalization by a factor of 5 to match the 
GRXE flux at the Carina Nebula. The GRXE emission 
contributes approximately the same flux as the CXB be- 
tween 1 and 8 keV, <10% of the total. 

Emission from LHB toward the Carina Nebula has 
an insignificant spatial variation over degrees, and has 
a surface brightness of ~4xl0 4 counts s _1 arcmin -2 in 
the ROSAT PSPC R1+R2 band (Snowden et al. 1998). 
Assuming a Raymond-Smith thin-thermal plasma model 
(Raymond & Smith 1977), the LHB emission with kT 
~0.1 keV accounts for ~30% of the 0.2-0.3 keV flux. 
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Fig. 5. Combined Chandra spectra of point sources in the north (left) and south (right) regions. The solid lines show the best-fit 
absorbed 2T models, with a Gaussian for the north spectrum. The dotted lines show each plasma component. The bottom panels 
show residuals from the best-fit models. 



Table 1. Best-fit Models of Chandra Point Source Spectra 



Parameters 




North 


South 




(keV) 


0.62 


0.27 


log E.M. cool * 


(cm -3 ) 


55.7 


55.9 


kT hot 


(keV) 


3.8 


4.2 


log E.M. hot * 


(cm -3 ) 


55.6 


55.6 


Abundance 


(solar) 


0.16 


0.08 




(10 21 cm -2 ) 


1.3 


1.7 


fW 


(10 -13 ergs cm -2 s -1 ) 


9.0 


7.5 



AxVd.o.f (d.o.f) 1.40 (63) 1.20 (51) 



* Assuming d =2.3 kpc, the distance to jj Carinae (Hillier et al. 2001). 

* Observed X-ray flux between 0.5-10 keV. 



Since emission from the Carina Nebula should be cut off 
around ~0.4 keV by interstellar absorption, the residual 
emission below 0.4 keV probably also arises in the fore- 
ground. 

Emission from the star i] Car is an order of magnitude 
stronger than any other source in the field and thus can 
contaminate a substantial fraction of the XIS FOV. We re- 
produced the rj Car spectrum using the same Suzaku data, 
modeling it with a commonly absorbed 2T model with an 
independent nitrogen abundance for the soft diffuse emis- 
sion below ~1 keV and a commonly absorbed 2T model 
for the central point source. We then estimated effective 
area of a source at the position of n Car on each source 
region using the arf generator xissimarf gen-2006-05-28 
(Figure 4) . The contamination is almost negligible for the 
south region and <10% for the north region. 

5.3. Total Contribution 

Figure 4 depicts the total contribution of potential con- 
taminants to the Suzaku spectra. Below 0.3 keV, the LHB 
accounts for ^30% of the observed flux. Above ~2 keV, 
the contamination sources contribute up to half of the 
Suzaku flux. This, in turn, means that half of the hard 
band flux cannot be attributed to the known contami- 



nants. On the other hand, the excess between 0.3 and 
2 keV is a factor of 3—10 higher than the total of possible 
contaminants and therefore seems to originate in diffuse 
plasma though a detailed analysis of the point source pop- 
ulation would be desirable. 

6. Spectral Fits 

To quantitatively estimate the spectral parameters of 
the diffuse emission, we simultaneously fitted the FI and 
BI spectra above 0.3 keV for each region. We tentatively 
assumed an absorbed IT thin-thermal plasma (WABS, 
Morrison & McCammon 1983; APEC code) for the HC, 
which should really originate from multiple emission com- 
ponents as seen in §5. For the MC, we first assumed an 
absorbed IT thin-thermal plasma, which however failed 
to reproduce both the spectral slope above 1 keV and 
the strong emission line from OVII ions. We therefore 
assumed a commonly absorbed 2T thin-thermal plasma 
(MC coo/ and MC hot ). We omitted data below 0.3 keV, 
that is the SC, and BI data points between 5.5—6.7 keV, 
where emission from 55 Fe calibration source dominates. 
We separately varied elemental abundances between the 
HC and MC. We also individually varied the abundance 
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Fig. 6. Best-fit result of the north (left) and south (right) spectra by the model with TVu(MC) and kT tied (black: Ff, red: Bf). 
The solid lines and dotted lines show the best-fit model and each of the plasma components, respectively. The bottom panels show 
residuals in the fits. 



of each element of the MC, but tied the abundances of the 
same elements between MC coo ; and MC hot- 

Calibration uncertainty is still a concern at the time of 
writing this paper. To adjust the effective area between 
the XIS sensors wc varied the BI normalization by multi- 
plying an energy independent coefficient to the Bl model 
("constant" model in xspec). To allow for uncertainty 
in the absolute energy scale of each spectrum we allowed 
the gain to vary ("gain fit" function in xspec). In the 
best-fit models, the BI spectra of the north and south 
regions have ~H% and ^5% higher normalization than 
the FI spectra, respectively. Other Suzaku observations 
also show ^3% normalization inconsistency between the 
BI and FI (e.g., NGC 2110, Bamba et al., private com- 
munication). A remaining ^5% discrepancy in the north 
spectra may be caused by inclusion of the 55 Fe calibration 
source region for XIS1 (~4% discrepancy in the extract- 
ing area). The slope and offset of the energy scale settled 
around ~2.1% and 12 eV for the BI south spectrum and 
^0.1—0.4% and ^1—7 eV for the others, values generally 
consistent with other rev 0.7 data. In subsequent fits, we 
fixed these parameters at their best-fit values. 

The best-fit models (Table 2) reproduce neither of the 
spectra with a goodness of fit acceptable above 90% con- 
fidence. This might be caused by poor calibration and/or 
inaccuracy in the emission line code, which we describe 
in detail later. In this model, plasma temperatures of 
both spectra are similar to each other, but the north 
spectrum has an order of magnitude larger cool compo- 
nent emission measure E.M.(MC coo i) and a factor of two 
larger Nn (MC) than the south spectrum. Fits of Chandra 
spectra also showed larger Ah to the north of the dif- 
fuse plasma (~7xl0 21 cm -2 ) than a negligible Ah to 
the south (Townsley 2006). However, factor of 2 larger 
Nh value for the north region seems inconsistent with 
a thick CO condensation being around the south region 
(Yonekura et al. 2005), though we do not know loca- 
tions of the X-ray plasma and the CO condensation along 
the line of sight. Moreover, it would be too convenient 
to think that the north region happens to have an or- 



der of magnitude larger E.M.(MC coo i), and, as a result, 
have a 0.3—0.5 keV observed flux similar to the south 
region. Since E.M.(MC CO oi) an d Ah(MC) can easily cou- 
ple in a spectral fit, it is possible that the differences of 
E.M.(MC coo i) and A H (MC) arise artificially as the result 
of a local x 2 minimum. 

We therefore assumed the north and south regions have 
the same absorption and tied together Ah(MC). We also 
tied together the temperatures kT(MG coo i), kT(MChot) 
and fcT^HC), based on their similarity in the uncon- 
strained fits. The best-fit model (Figure 6, Table 2) is once 
again not acceptable above 90% confidence, and the resid- 
uals are similar to these in the individual fits. A broad 
residual hump is visible around 0.65—0.80 keV, where 
emission from the Fe-L shell complex is strong. The hump 
is not reduced using a multi-temperature model (c6pvmkl) 
in xspec. This may mean that the plasma is not in col- 
lisional equilibrium, or that the emission line code does 
not reproduce the Fc L-shell line complex. On the other 
hand, the model does not reproduce the argon and sulfur 
lines around 2—3 keV, confirming that the HC has mul- 
tiple components (see §5). Both spectra have marginal 
peaks with equivalent width of ~50— 60 eV at 6.4 keV, 
which might arise from fluorescence of cold iron. 

Although the best-fit models are not formally accept- 
able, they seem to reproduce the global spectral shapes 
well. The plasma temperatures, ~0.2 and ~0.6 keV 
in the MC and ~5.5 keV in the HC, arc similar to 
the results of the individual fits. The E.M.(MC coo i) 
and E.M.(MChot) have the same ratio in the north and 
south spectra, confirming that the two spectra are quite 
similar between 0.3 and 2 keV. The average surface 
fluxes of the north and south spectra are 7 and 9.5 
xlO -14 ergs cm -2 s _1 arcmin -2 between 0.5—2 keV and 
3.5 and 3.8 xlO -14 ergs cm -2 s" 1 arcmin~ 2 between 
2—8 keV, respectively. The 0.5—2 keV band surface flux 
is consistent with the unabsorbed surface flux estimated 
from the least contaminated Chandra spectrum, assum- 
ing N R =3xl0 21 cm~ 2 (Evans et al. 2003). This result 
ensures that the Suzaku spectra between 0.3—2 keV rep- 
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Model 



Individual 
North South 



ATh(MC), kT tied 
North South 



JVh(MC) fixed NEI Typical Error* 

North South North South 



Medium Component (MC) 






















(keV) 




0.20 


0.23 




0.21 


0.19 




0.21 


0.01 


kT hot 


(keV) 




0.60 


0.58 




0.58 


0.56 




0.57 


0.01 


log E.M. coo J 


(cm _3 arcmin~ 


2 ) 


56.0 


55.0 


55.7 


55.6 


56.2 


56.0 


55.5 55.3 


0.1 


log E.M. ho J 


(cm _3 arcmin _ 


2 ) 


55.0 


54.9 


55.1 


55.0 


55.0 


54.8 


55.1 55.0 


0.1 


N H 


(10 21 cm- 2 ) 




2.4 


1.2 




1.8 


3.0 (fix) 


1.8 


0.2 


C 


(solar) 




2.5e-2 


0.0 


0.0 


0.0 


0.34 


0.85 


0.0 0.25 


5e-2 


N 


(solar) 




0.0 


1.2e-2 


0.0 


0.0 


0.0 


0.0 


1.5e-2 4.5e-2 


lc-2 





(solar) 




4.3e-2 


0.12 


6.2e-2 


6.3e-2 


4.2e-2 6.1e-2 


9.1e-2 0.12 


le-2 


Ne 


(solar) 




8.9e-2 


0.14 


0.13 


9.0e-2 


8.1e-2 


0.11 


0.19 0.14 


3e-2 


Mg 


(solar) 




0.11 


0.19 


0.11 


0.16 


9.4c-2 


0.17 


0.11 0.16 


2e-2 


Al 


(solar) 




0.17 


0.0 


7.5e-2 


2.9e-2 


0.17 


0.18 




0.2 


Si 


(solar) 




0.10 


0.25 


9.2e-2 


0.24 


0.11 


0.36 


0.11 0.27 


2e-2 


S 


(solar) 




0.38 


0.12 


0.24 


0.39 


0.48 


0.85 


0.22 0.43 


0.1 


Ar 


(solar) 




0.36 


0.0 


0.0 


0.13 


1.4e-2 


1.2 


1.4e-2 1.2 


0.2 


Ca 


(solar) 




0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 0.0 


3e-2 


Fe 


(solar) 




8.4c-2 


0.30 


7.8c-2 


0.30 


8.3e-2 


0.43 


8.3e-2 0.30 


le-2 


Ni 


(solar) 




0.34 


0.56 


8.9e-2 


0.78 


0.56 


1.9 


1.5e-2 0.39 


0.1 


T 


(s cm -3 ) 
















1.7cl2 


> 9cll 


Hard Component (HC) 




















kT 


(keV) 




5.5 


5.5 




5.3 


5.5 




5.1 


0.2 


log E.MJ 


(cm- 3 arcmin- 


2 ) 


54.4 


54.5 


54.4 


54.5 


54.4 


54.5 


54.5 54.5 


0.1 


Nh 


(10 22 cm" 2 ) 




2.3 


1.3 


2.4 


1.5 


1.8 


1.2 


2.5 1.7 


1 


Z 


(solar) 




0.48 


0.44 


0.47 


0.43 


0.50 


0.44 


0.46 0.42 


5e-2 


X'fd.o.f 






2.07 


2.94 




2.70 


3.52 




2.70 




d.o.f. 






378 


344 




736 


737 




737 





Assuming d =2.3 kpc. 



A X 2 =2.7. 



resent diffuse emission from the Carina Nebula. 

The south region has a factor of 2 higher Si abundance 
and a factor of 4 higher iron abundance than the north 
region. Abundances of magnesium and nickel may also be 
enhanced by a factor of ~1.5 and 7—8 in the south, though 
these lines are situated among the iron L-shell lines and 
their abundance measurements are strongly affected by 
the remaining uncertainty in fitting the iron line complex. 
No other elements show a significant abundance difference. 

The derived Nu (^1.8xl0 21 cm -2 ) is consistent with 
ATh to some objects around the Carina Nebula (Seward 
et al. 1979; Savage et al. 1977). However, Leutenegger, 
Kahn and Ramsay (2003) argued that these measurements 
are problematic and suggested a factor of 1.5 higher N-r 
(~3xl0 21 cm- 2 ). When we fixed N H at 3xl0 21 cm- 2 in 
the spectral fit, the abundances of silicon, sulfur and car- 
bon, which have emission lines in the 0.4—0.5 keV band 
where the spectrum cuts off sharply due to interstellar ab- 
sorption, significantly increase (see Table 2). This means 
that abundances of these elements strongly couple with 
./Vh, and are therefore less reliable. 

Since some emission lines such as MgXI and FeXVII 
have residuals on the low energy side, we also tried a non- 
equilibrium ionization (NEI) model (Table 2). However, 



the large ionization time scale r ^1.7xl0 12 s cm -3 in the 
best-fit rather suggests the ionization equilibrium plasma, 
and therefore the best-fit model did not reduce the resid- 
uals. We also tried multi-temperature models for the cool 
component, but the result did not improve, either. 

7. Discussion 

Based on the morphology of the X-ray emission, the X- 
ray luminosity (~10 35 ergs s" 1 ) and the plasma tempera- 
ture (~0.8 keV), Seward & Chlebowski (1982) argued that 
diffuse X-ray emission from the Carina Nebula is proba- 
bly produced by strong stellar winds from early stars col- 
liding with ambient gas. Our results using Suzaku and 
XMM-Newton are basically consistent with these previ- 
ous results, except for a slightly lower temperature (~0.2, 
0.6 keV). The inferred temperature is also similar to the 
temperature of the diffuse emission of the Omega and 
Rosette Nebulae, for which Townsley et al. (2003) sug- 
gested a stellar wind origin. 

What is new in this result is the detailed line diagnos- 
tics in the soft energy band, especially below ~1 keV. 
The study reveals that nitrogen is significantly depleted 
relative to other heavy elements throughout the observed 
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regions, while the abundance of iron and silicon in the 
southern portion of the nebula is 2-4 times higher than in 
the northern portion. Any successful explanation of the 
diffuse emission from the Carina Nebula must account for 
these abundance anomalies. In the paragraphs below, we 
explore possible mechanisms. 

The evolved stars in the field, rj Car and WR 25, have 
strong winds rich in nitrogen and poor in oxygen, result- 
ing from the CNO hydrogen burning process (Davidson 
et al. 1982; van der Hucht et al. 1981). Strong nitro- 
gen and negligible oxygen lines are observed from the 
diffuse X-ray plasma immediately surrounding r\ Car ; 
the abundance ratio was restricted to N/O >9 (Tsuboi 
et al. 1998; Leutenegger et al. 2003; and the -q Car spec- 
trum from the same Suzaku observation). Although less 
evolved massive O stars in general show a nitrogen over- 
abundance, the O stars in the Carina region have a rela- 
tively low N/O ratio (Morrell et al. 2005). Even so, the 
Suzaku spectra of the Carina Nebula show no significant 
nitrogen emission, but show clear oxygen emission lines. 
The N/O abundance ratio inferred from the spectral fits 
is <0.4, over 20 times less than around r\ Car. The abun- 
dance distribution is totally contrary to that expected 
from massive stellar winds, unless the winds somehow heat 
the interstellar matter without enriching it, thus leaving 
the X-ray plasma with abundances typical of interstellar 
matter. At the same time, X-ray luminosity of the Carina 
Nebula is about two orders of magnitude higher than that 
of other Galactic star forming regions, but the number of 
early O stars is only an order of magnitude higher (see 
Table 4 in Townsley et al. 2003). These results suggest 
an additional energy source is needed to power the X-ray 
emission in the Carina Nebula. 

An obvious possibility is one or more core-collapse su- 
pcrnovae (i.e. Type Ib,c or II), mentioned as a possi- 
bility by Townsley et al. (2003). The north and south 
regions vary strongly in silicon and iron abundances 
and marginally in magnesium and nickel abundances. 
These elements are products of core-collapse supernovae. 
Moreover, young SNRs such as Cas A and Vela show 
strong abundance variation from location to location (e.g., 
Katsuda & Tsunemi 2006). 

A single supernova could power the entire Carina 
Nebula. Assuming a uniform temperature throughout the 
entire extent of the nebula (l.°l E-W by 0.°7 N-S (Seward 
& Chlebowski 1982) by a comparable line-of-sight depth) 
the total energy content in the hot gas of ~2xl0 50 ergs, 
a modest fraction of the ^10 51 ergs of kinetic energy pro- 
duced by a canonical supernova. 

On the other hand, the total iron mass in the diffuse gas 
cannot be the product of a single supernova. Assuming an 
iron abundance of 0.30 solar throughout the entire nebula 
implies a total iron mass Mp e ~ 0.1/n M Q , where n is the 
plasma number density in the unit of cm~ 3 . Our results 
suggest n ~0.2— 0.4 cm -3 in the north and south regions, 
assuming the diffuse plasma in each region has a dimen- 
sion along the line of sight comparable to its project angu- 
lar extent at a distance of 2.3 kpc. If these regions are rep- 
resentative of the entire nebula then Mp e ~0.25— 0.4 Mq. 



The typical core collapse supernova yields <0.1 M Q , so 
more at least 3-5 supernovae would be required to supply 
the iron. 

There is no evidence for supernova remnants in the 
Carina Nebula in the radio (Whiteoak 1994) or the X-ray. 
There are two known pulsars, IE 1048.1-5937 and PSR 
J1052— 5954, within 1° from r] Car, but they are located 
far behind the Carina Nebula (Kramer et al. 2003; Gavriil 
& Kaspi 2004; Gaensler et al. 2005). Core collapse super- 
novae produce a wide variety of stellar remnants, however, 
and not all of them pulse. It is plausible that several stel- 
lar remnants similar to those in Cas A and Puppis A (with 
Lx ~2xl0 33 ergs s -1 and fcT~5xl0 6 K) could reside un- 
detected in the Carina Nebula. There is no requirement 
that the stellar remnants are located within the north or 
south regions for which we carefully inspected the point 
sources, or even in the Carina Nebula after ^10 6 yr. 

The need for multiple supernovae to explain the abun- 
dances with the modest energy input requirements is 
consistent with a hypothesis advanced by Yonekura 
et al. (2005) to account for the fact that the strong turbu- 
lence in the r/ Car giant molecular cloud (GMC) observed 
in CO lines is too energetic to be produced by stellar winds 
or a single SNR. They suggest that the turbulence was 
preexisting when the GMC was formed. They postulate 
that it was produced by the large number of supernovae 
(>20) responsible for creating the Carina flare supershcll 
(Fukui et al. 1999), and is now dissipating. On the other 
hand, Smith et al. (2000) found a bipolar morphology in 
the Carina Nebula in the infrared and suggested the pres- 
ence of one or more supernova remnants. Whichever case 
is connected to the diffuse X-ray emission, the absence of 
evidence of SNRs in the Carina Nebula is not a problem, 
as the explosions that power the X-ray nebula (and the 
molecular cloud turbulence) occurred ~10 6 years ago. 

The abundance variation between the north and south 
regions could be one vestige of the multiple explosions. 
Alternatively, since silicon, iron and magnesium are re- 
fractory elements required for silicate grains and can be 
very strongly depleted from the gas phase in the ISM 
(Salpetcr 1977; Mathis 1996; Draine 2003), we may con- 
sider the dust formation or destruction, both of which are 
strongly associated with the environment of evolved mas- 
sive stars, to explain the abundance variation. 

SNRs thus seem to be a promising hypothesis for 
explaining the diffuse plasma from the Carina Nebula. 
However, there are several issues yet to be addressed to 
confirm this hypothesis, (i). The absolute abundances of 
all the heavy elements are smaller than solar, unlike in a 
canonical supernova remnant. Other Suzaku observations 
of diffuse plasma, including the Cygnus Loop supernova 
remnant (Miyata et al. 1998), also show similar low ab- 
solute elemental abundances. This could be simply due 
to a problem of the current X-ray measurement, such as 
the need to include a non-thermal continuum component, 
unknown physical processes of line emission in a very low 
density plasma, or helium-rich plasma, (ii). The oxygen 
abundance is relatively low though it should be the most 
abundant element in core collapse supernova ejecta. In 
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particular, the O/Fe ratio is smaller than the solar ratio. 

What is the excess below the OVII line in the north 
spectrum? Since the other emission lines did not show 
such an apparent broadening, poor energy calibration is 
an unlikely cause. The excess can be reproduced by a nar- 
row Gaussian line with a center energy of 0.547±0.01 keV. 
The error range includes three emission lines from SXIV 
(10->1, 11->1) and CaXVI (95->8). However, sulfur 
and calcium should have stronger peaks below <0.5 keV 
when kT ~0. 1—0.8 keV, which are not seen in the spec- 
trum. Though strong radiation from Tr 14 cluster and 
rj Car could photo-ionize the plasma, the excess peak is 
slightly lower than the forbidden line of O VII at 0.561 keV, 
which is strong in the photo-ionized plasma (Porquet & 
Dubau 2000). If the excess is produced by the Doppler 
shift of the helium-like oxygen line, the velocity is ex- 
tremely high (^7.5xl0 3 km s _1 ). We have no good ex- 
planation for this feature. We should note that a similar 
excess below O VII is seen from a few Suzaku observations 
such as the North Polar Spur (Miller et al. in prepara- 
tion) . 

Less than half of the surface flux from the hard com- 
ponent arises from known contaminants: X-ray point 
sources, the CXB and the GRXE. One possible ori- 
gin of the remainder is a large number of low mass 
young stars embedded in the cloud, which are too 
faint to be detected individually. Both the north 
and south regions include massive clouds (Yonekura 
ct al. 2005), which probably have the star forming 
activity. Low mass pre-main-sequence stars in the 
Orion Nebula emit X-rays of fcT~0.5— 3 keV and Lx 
_ 10 28.9 ergs s -i between 2-8 keV (Feigelson et al. 2005). 
Since 2—8 keV X-ray surface brightness of the Carina 
Nebula is 2.3 xlO 31 ergs s _1 arcmin -2 , the surface num- 
ber density of young stars should be ~290 arcmin -2 = 
650 pc -2 at 2.3 kpc. Since the star formation would occur 
over a parsec scale along the line of sight, the number den- 
sity should be <650 pc -3 , which is small compared with 
the Orion Nebula (peak of >10 4 pc" 3 , Hillenbrand 1997). 
Multiple X-ray sources detected in a deep Chandra obser- 
vation of the Tr 14 cluster (Townsley 2006) may account 
for a part (or all) of this possible hidden population. There 
are two problems: hard X-ray surface flux is almost the 
same in the two regions even though there are more po- 
tential star formation sites in the north; and kT ^5.5 keV 
is higher than the typical temperature of low mass young 
stars though a possibly different population in the Carina 
Nebula might produce a composite spectrum with a higher 
effective plasma temperature. Another possibility is that 
the hard X-ray emission comes from diffuse plasma. Such 
a very high temperature plasma perhaps needs young su- 
pernova remnants (e.g., Koyama et al. 2006a), though the 
observed soft X-ray emission rather suggests presence of 
old supernova remnants. 

8. Summary 

We analyzed diffuse X-ray emission from the Carina 
Nebula using the CCD camera XIS onboard the Suzaku 



observatory The XIS has the best soft band spectral reso- 
lution for extended sources of any X-ray observatory, and 
provides detailed line profiles and reliable elemental abun- 
dances of the plasma. To understand spatial distribution 
of the diffuse emission, we first analyzed an XMM-Newton 
MOS image centered at rj Car. An X-ray color map in- 
dicates that the diffuse emission north of rj Car is softer 
than the emission to the south. Analysis of the Suzaku 
spectra shows that this difference arises from a differ- 
ence in the iron L-shcll line intensity. Suzaku spectra are 
roughly divided into three X-ray components. Spectra 
below ~0.3 keV have almost the same surface flux in the 
north and south regions. This component originates in the 
foreground, including the LHB. The spectra above 2 keV 
also have almost the same surface flux and shape. About 
half of the emission arises from the CXB, GRXE and X- 
ray point sources detected with Chandra. The remainder 
of the hard emission may be from unresolved pre-main- 
sequence stars. The spectra between 0.3 and 2 keV are 
dominated by emission from soft diffuse emission. The 
spectra show K-shell lines of silicon, magnesium, neon, 
oxygen, and possibly carbon ions, and L-shell lines of iron 
ions. The strength of the Si and Fe lines is different be- 
tween the north and south regions, while the other parts 
of the spectra are similar. The spectra from both regions 
are reproduced with a 2T plasma model with fcT^0.2 and 
0.6 keV and sub-solar elemental abundances and common 
absorption of ~2xl0 21 cm -2 , but are not acceptable at 
the 90% confidence level. No nitrogen line is detected from 
either of the spectra, restricting N/O <0.4. Such a low ni- 
trogen abundance is not expected from stellar winds from 
evolved massive stars, and suggests that the plasma origi- 
nates from an old SNR, or a super shell produced by mul- 
tiple SNRs. The south region has a factor of 2—4 higher 
elemental abundances in iron and silicon than the north 
region. Magnesium and nickel may show spatial abun- 
dance variation, as well. The abundance variation may be 
explained by the presence of SNR ejecta. Alternatively, 
dust formed around the star forming core may play an 
important role. 

The spatial structure and abundance variation of the 
plasma provide clues to the origin of the diffuse emis- 
sion. For this purpose, a mapping of the entire Carina 
Nebula with Suzaku and XMM-Newton will be important. 
Meanwhile, the population of faint stars in the Carina 
Nebula should be investigated with a deep Chandra obser- 
vation, through a study of the X-ray luminosity function. 
The spectra are complicated with many emission lines; 
high resolution spectroscopy on the future satellites, such 
as an X-ray calorimeter, will be needed to solve the origin 
of this mysterious X-ray nebula. 
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